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ABSTRACT 

Context. The detailed chemical abundances of extremely metal-poor (EMP) stars are key guides to understanding the early chem- 
ical evolution of the Galaxy. Most existing data are, however, for giant stars which may have experienced internal mixing later. 
Aims. We aim to compare the results for giants with new, accurate abundances for all observable elements in 18 EMP turnoff stars. 
Methods. VLT/UVES spectra a.t R ~ 45,000 and S/N~ 130 per pixel {AA 330-1000 nm) are analysed with OSMARCS model 
atmospheres and the TURBOSPECTRUM code to derive abundances for C, Mg, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, Zn, Sr, and Ba. 
Results. For Ca, Ni, Sr, and Ba, we find excellent consistency with our earlier sample of EMP giants, at all metallicities. However, 
our abundances of C, Sc, Ti, Cr, Mn and Co are ~0.2 dex larger than in giants of similar metallicity. Mg and Si abundances are 
~0.2 dex lower ( the giant [Mg/Fe] values are slightly revised), while Zn is again ~0.4 dex higher than in giants of similar [Fe/H] 
(6 stars only). 

Conclusions. For C, the dwarf/giant discrepancy could possibly have an astrophysical cause, but for the other elements it must 
arise from shortcomings in the analysis. Approximate computations of granulation (3D) effects yield smaller corrections for giants 
than for dwarfs, but suggest that this is an unlikely explanation, except perhaps for C, Cr, and Mn. NLTE computations for Na 
and Al provide consistent abundances between dwarfs and giants, unlike the LTE results, and would be highly desirable for the 
other discrepant elements as well. Meanwhile, we recommend using the giant abundances as reference data for Galactic chemical 
evolution models. 

Key words. Galaxy: abundances - Galaxy: halo - Galaxy: evolution - Stars: abundances - Stars: Population II - Stars: Supernovae 



* Based on observations obtained with the ESO Very Large Stars", ID 165.N-0276; PL: R. Cayrel, and Programme 078.B-0238; 
Telescope at Paranal Observatory, Chile (Large Programme "First PL: M. Spite). 
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1. Introduction 

The surface composition of a cool star is a good diagnostic of 
the chemical composition of the gas from which it formed, if 
mixing with material processed inside the star itself has not oc- 
curred. Cool, long-lived stars have thus been extensively used 
to study the early chemical evolution of our Galaxy (and, by 
impUcation, other galaxies as well). The trends in abundance 
ratios which have been established over the last 30 years pro- 
vide important constraints on the early chemical evolution of 
the MiUcy Way (see Cayrel 1 1 9961 120061 for classic and recent 
reviews of the topic). 

Our own programme, "First Stars", is a comprehensive 
spectroscopic study of extremely metal-poor (BMP) stars to 
obtain precise information on the chemical composition of the 
early ISM and the yields of the first generation(s) of super- 
novae, conducted with the VLT and UVES spectrograph. The 
target stars have been selected from the medium-resolution 
follow-up (Beers et al. in preparation; Allende Prieto et al. 
I2000J of the HK objective-prism survey (Beers et al. 119851 
1 19921 and Beers [T999] l. initiated by George Preston and Steve 
Shectman, and later substantially extended and followed up by 
Beers as part of many collaborations, including the present one. 

Several papers have presented our results on the giant stars, 
which lend themselves most readily to the study of many ele- 
ments: Hill et al. ( |2002| - First Stars I), Depagne et al. (|^02| - 
First Stars II), Frangois et al. (12003] - F irst St ars III), Cayrel et 
al. (2004 - Firs t Stars V), Spite et al. ( [20051 - First Stars VI), 
Frangois et al. dlOOTl - First Stars VIII), and Spite et al. (2006 
- First Stars IX). In these papers, we found the abundances in 
some giants to have been altered with respect to their initial 
chemical composition, due to mixing with layers affected by 
nuclear burning. All the stars have undergone the first dredge- 
up, so their abundances of Li, C, and N are under suspicion. 
However, our detailed analysis (First Stars VI and IX), showed 
that the surface abundances of the less luminous giants (those 
below the "bump" in the luminosity function) are not signifi- 
cantly affected by mixing. 

It is therefore expected that the less-luminous giants and 
dwarfs should display the same abundances, provided that the 
surface composition of the latter has not been changed by atmo- 
spheric phenomena, such as diffusion. Comparing abundance 
ratios in dwarfs and giants can therefore, in principle, yield in- 
sight into the degree of mixing in giants and diffusion in dwarfs 
as well as which element ratios are reliable guides to the com- 
position of the early ISM in the Galaxy. 

So far, only few of our papers have discussed results for 
BMP dwarfs: Sivarani et al. (|2U04l - Fii'st Stars IV. .2006, - First 
Stars X), Bonifaci o et a l. (120071 - First Stars VII), and Gonzalez 
Hernandez et al. (2008'- First Stars XI). First Stars VII focused 
on the Li abundance, but also discussed the model parameters 
and [Fe/H] of the dwarf sample in considerable detail. Here 
we discuss the abundances from C to Ba in the same stars and 
compare the results for dwarfs and giants. 

2. Observations and reduction 

The sample of stars and the observational data are the same 
as discussed in Paper VII (Bonifacio et al. 120071 ). The ob- 



servations were performed with the ESO VLT and the high- 
resolution spectrograph UVES (Dekker et al. 120001 ) at a reso- 
lution of R- 45,000 and typical S/N ratios per pixel of -130 
on the coadded spectra (average 5 pixels per resolution el- 
ement). The spectra were reduced using the UVES context 
within MIDAS (Ballester et al. 120001) : see paper V for details. 
The region of the Mg i b triplet in our spectra is shown in Fig. 
[T](see also Fig. 1 of Paper VII, which shows the Li line in the 
same stars). Equivalent widths were measured on the coadded 
spectra. For a few stars, for which spectra with different resolu- 
tions (different slit-width, or image sheer used) were available, 
we coadded separately the spectra with the same resolution and 
then averaged the equivalent widths. 

3. Determination of atmospheric parameters 

We have carried out a classical ID LTE analysis using 
OSMARCS models (see, e.g., Gustafsson et al. [T9751 120031 
12008 ). Estimates of Tgff were derived from the wings of Ha; 
log g estimates were obtained by consideration of the ionisa- 
tion equilibrium of iron. Microturbulent velocities were fixed 
by requiring no trend of [Fe i/H] with equivalent width. Details 
are given in "First Stars VII", together with an extensive dis- 
cussion of the effective temperature scale. In that paper we es- 
tablished that our Ha based temperatures satisfy the iron ex- 
citation equilibrium and are also in good agreement with the 
Alonso et al. ( 1996 ) colour-temperature calibration, which we 
used for the giant stars (Cayrel et al. 120041 ) The adopted param- 
eters are listed in Table [1] 

The parameters of the subgiant star BS 16076-006 require 
a comment, because the Balmer line broadening in this star in- 
creases from Ha towards the higher members of the Balmer 
series. Our adopted Ten (5199K) is derived from the wings of 
Ha, but the wings of HS correspond to a much higher effective 
temperature, of the order of 5900 K. All values of Tgff derived 
from colours are also consistently higher than derived from the 
Ha profile, confirming this peculiarity. This star was also anal- 
ysed from medium-resolution ESI - Keck spectra (R-IQQQ) by 
Lai et al. (2004), who adopted a TefF= 5458 K, based on pho- 
tometry. Such a TefF is compatible with the profile of Hy, but too 
low to reproduce the profile of HS. The reason for this peculiar 
behaviour (e.g. a binary companion or chromospheric activity), 
needs further investigation, but the three radial velocities of de- 
rived from our two spectra and that of Lai et al. (2004) show 
no evidence of variation. None of our results depends critically 
on the abundances of this star, however. 



4. Abundance determination 

The abundance analysis was performed using the LTE spec- 
tral line analysis code "Turbospectrum" (Alvarez and Plez, 
1998 ). The abundances of the different elements have been de- 
termined mainly from the equivalent widths of unblended lines. 
However, synthetic spectra have been used to determine abun- 
dances from the molecular bands, or in cases when the lines 
were severely blended, affected by hyperfine structure, or were 
strong enough to show significant damping wings (see Sect. 
I6.lt . The abundances of C and the a elements (as well as for 
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Table 1. Adopted model atmosphere parameters. Our UVES 
spectra show that BS 16076-006 is in fact a subgiant (First Stars 
VII), and CS 29527-015 a double-lined binary; these stars are 
omitted in Fig.[T] All the others seem to be single turnoff stars. 
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Sc) are listed in Table |2l those of the heavier (neutron-capture) 
elements are listed in Table [3] 

Abundance uncertainties are discussed in detail in Cayrel 
et al. (2004) and Bonifacio et al. ( |2007| l. For a given temper- 
ature, the ionisation equilibrium provides an estimate of the 
gravity with an internal precision of about 0. 1 dex in log g, and 
the microturbulent velocity can be constrained within about 0.2 
km s ' . The largest uncertainty comes from the temperature de- 
termination, which is uncertain by ~100K. 

The total error estimate is not the quadratic sum of the var- 
ious sources of uncertainty, because the covariance terms are 
important. As an illustration of the total expected uncertainty 
we have computed the abundances of CS 29177-009 with dif- 
ferent models: Model A has the nominal temperature 6260 K, 
gravity (log g - 4.5), and microturbulent velocity (v, = 1.3 
km s '), while Models B and C differ in log g and v, by Icr. 
Model D has a temperature 100 K lower and the same log g 
and V,, while in Model E we have determined the "best" values 
of log g and v, corresponding to the lower temperature. The 
detailed results of these computations are given in Table |4] 

5. C, N, O abundances 

5.7. Carbon 

The carbon abundance was determined by spectrum synthesis 
of the A^A - X^n band of CH (the G band ). Wav elengths of 
the CH lines are from Luque and Crosley ( 119991 ); transition 
energies are from the list of J0rgensen et al. (1996); isotopic 
shifts were computed using the best set of available molecular 
constants. The strongest lines of '^CH at 423nm are invisible in 
all of our stars, so the '^C/'^C ratio could not be measured. In 
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Fig. 1. The region of the Mg i b triplet in the program stars. 
[Fe/H] is shown to the left of each spectrum. In these EMP 
stars, the triplet lines have no damping wings. 
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Fig. 2. [C/Fe] ratios in our turnoff stars (black dots) and un- 
mixed giants (open circles). The grey triangle shows the sub- 
giant BS 16076-006. 




Fig. 3. [N/Fe] in our sample of unmixed giants. The triangle at 
[Fe/H]= -3.81 shows the subgiant BS 16076-006. 



computing the total C abundance, we have therefore assumed a 
solai- '^c/i^c ratio. 

In Fig. |2] we present the measured [C/Fe] values in our 
dwarf stars and compare them to values for our unmixed giants 
from Paper V. In this figure we have omitted the mixed giants, 
located above the bump, since we have shown (First Stars VI 
and IX) that the abundances of C and N in the atmospheres 
of these stars are strongly affected by mixing and thus are not 
good diagnostics of their initial chemical compositions. 

The mean [C/Fe] value for the turnoff stars is [C/Fe] - 
0.45 + 0.10 (s.d.), but [C/Fe] = 0.19 + 0.16 (s.d.) for the giants. 
Thus, we find a moderately significant difference between the 
C abundances in the giants and the turnoff stars (Fig|2]l. We 
discuss the possible causes of this discrepancy in SecHT] The 
mean [C/Fe] has been computed excluding the binary turnoff 
star CS 29527-015, which appears to be quite carbon rich (Fig. 



Hernandez et al. 120081 Paper XII), and these are compatible 
with the O abundances measured in giants. 

Our spectra of the dwarfs discussed in the present paper do 
not cover the OH lines in the UV. The success in the case of 
CS 22876-032 suggests that these lines probably offer the only 
option for measuring O abundances in EMP dwarfs. 

6. The a elements: Mg, Si, Ca, Ti 

Fig. ID presents the observed [a/Fe] ratios in our EMP dwarf 
and giant samples. 

A priori, we expect to find the same mean abundance for 
these elements in dwarfs and in giants, and this is what we see 
for Ca. However, the mean [Mg/Fe] and [Si/Fe] ratios are ~ 0.2 
and 0.3 dex lower in the EMP dwarfs than in the giants, while 
the mean abundance of [Ti/Fe] ratio is higher in the dwarfs 
by about 0.2 dex. What are the possible causes of these differ- 
ences? 



5.2. Nitrogen 

Generally, the NH (and CN) bands are not visible in the spectra 
of EMP turnoff stars (the stars are too hot), so N abundances 
can only be measured in strongly N-enhanced stars (First Stars 
X). The subgiant BS 16076-006 exhibits a weak NH band, 
however, and we find [N/Fe]= +0.29 for this star, taking into 
account the correction of -0.4 dex derived in Paper VI. 

Figure [3] shows the measured [N/Fe] ratios for our sam- 
ple of "unmixed" giants (Paper VI). BS 16076-006 agrees with 
(and thus supports) the high [N/Fe] values found in the giants 
at the lowest metallicities. 



5.3. Oxygen 

We have not been able to measure O abundances for any of our 
dwarf stars. The [OI] Une at 630.03 nm, which we used for 
giants, is too weak, as is the permitted O i triplet at 770 nm, 
given the S/N we achieve in this spectral region. Only for the 
dwarf binary system CS 22876-032 have we been able to mea- 
sure O abundances, using the OH lines in the UV (Gonzalez 



6.1. Magnesium 

In Fig. in the Mg abundance for the giant stars has been derived 
from a full fit to the profiles of the Mg lines, in contrast to the 
results given by Cayrel et al. j2004l Paper V). The equivalent 
widths of the Mg lines are often quite large (EW > 120mA), 
and in Paper V we underestimated the equivalent widths of 
these lines by neglecting the wings. For the most Mg-poor stars 
in our sample the lines are weak and the difference negligible, 
but it is quite significant in most of our stars, with a mean sys- 
tematic difference of about 0.15 dex. In the dwarfs, the abun- 
dance has been derived from profile fits to the strongest lines 
(the lines at ~383 nm, which are also located in the wings of a 
Balmer line); and from equivalent widths for the weak lines. 

6.2. Silicon 

In the cool giants the Si abundance is derived from a line at 
410.3 nm. This line (multiplet 2) is located in the wing of Hd, 
and the hydrogen line has been included in the computations. 
There is another line at 390.6 nm (multiplet 3), but in giants this 
line is severely blended by CH lines. In turnoff stars the line at 
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Table 3. Abundance ratios for the iron-peak and neutron-capture elements. 
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+0.31 


1 






9 


CS 22966-011 


-3.07 


-0.23 


0.03 5 


-0.70 


0.00 3 


0.48 


0.12 


4 


0.08 


0.06 2 


< 0.50 




0.03 


1 


-0.05 


1 


10 


CS 29499-060 


-2.70 


0.01 


0.04 6 


-0.28 


0.02 3 


0.36 


0.09 


4 


0.19 


0.09 3 


0.73 


1 


-0.60 


1 






11 


CS 29506-007 


-2.91 


-0.12 


0.05 5 


-0.59 


0.01 3 


0.39 


0.03 


3 


0.04 


0.08 3 


0.71 


1 


0.16 


1 


0.18 


1 


12 


CS 29506-090 


-2.83 


-0.16 


0.06 5 


-0.62 


0.02 3 


0.45 


0.11 


4 


0.04 


0.12 3 


0.66 


1 


0.36 


1 


-0.35 


1 


13 


CS 29518-020 


-2.77 


-0.18 


0.05 2 












0.04 


- 1 


<0.33 












14 


CS 29518-043 


-3.20 


-0.20 


0.08 4 


-0.64 


0.00 2 


0.57 




1 


0.07 


0.01 2 


<0.68 




0.08 


1 






15 


CS 29527-015 


-3.55 


-0.21 


0.15 4 


-0.66 


- 1 


0.70 




1 


-0.09 


0.05 2 


<0.98 




0.34 


1 






16 


CS 30301-024 


-2.75 


-0.16 


0.06 5 


-0.59 


0.01 3 


0.30 


0.11 


4 


0.02 


0.04 3 


0.55 


1 


-0.32 


1 


-0.28 


1 


17 


CS 30339-069 


-3.08 


-0.24 


0.06 5 


-0.71 


0.00 3 


0.33 


0.05 


2 


-0.01 


0.17 3 


< 0.47 




-0.10 


1 






18 


CS 31061-032 


-2.58 


-0.10 


0.16 6 


-0.51 


0.02 3 


0.38 


0.15 


4 


0.03 


0.05 3 


0.40 


1 


0.21 


1 


-0.40 


1 




BS 16076-006 


-3.81 


-0.41 


0.16 6 


-0.93 


0.10 3 


0.39 


0.05 4 


-0.05 


0.04 3 






<-1.59 


1 


<-1.0 


1 



[Mn/Fe] has been determined only from the lines of the resonance triplet. 



410.3 nm is invisible, but the CH lines are weak enough that 
the line at 390.6 nm can be used. Thus, in the end, only a sin- 
gle Si line (but not the same one) could be used in both dwarfs 
and giants; a systematic error in the log of these lines could 
explain the observed difference. Both lines are in fact measured 
in the subgiant star BS 16076-006 and yield consistent Si abun- 
dances, but, given the uncertain atmospheric parameters of this 
star (see[3]l, a systematic error in log cannot be ruled out. Our 
new [Si/Fe] ratios are in good agreement with the value found 



from the same Si line by Cohen (I2004I I. also for EMP turnofF 
stars. 



6.3. Titanium 

The Ti I lines are very weak in turnofF stars, so the Ti abun- 
dance can only be determined from the Ti II lines. About 15 
Ti II lines could be used, and the internal error of the mean 
is very small (less than 0.1 dex). Fig. H] clearly shows higher 
[Ti/Fe] ratios in the dwarfs than in the giants (A[Ti/Fe] = 0.2 
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Table 4. Abundance uncertainties linked to stellar parameters. 



CS 22177-009 

A: Teff=6260 K, log g = 4.5, vt=1.3 km s" 
B: Teff=6260 K, log g = 4.4, vt=1.3 km s" 
C: Teff=6260 K, log g = 4.5, vt=l.l km s" 
D: Teff=6I60 K, log g = 4.5, vt=1.3 km s" 
E: Teff=6160 K, log g = 4.3, vt=1.2 km s" 



El. 




Ac-A 


Ad-a 


Ae-a 


[Fe/H] 


-0.01 


0.03 


-0.05 


-0.06 


[Na I/Fe] 


0.02 


-0.02 


-0.01 


0.01 


[Mg I/Fe] 


0.03 


-0.01 


-0.02 


0.00 


[Al I/Fe] 


0.01 


-0.03 


-0.03 


-0.01 


[Si I/Fe] 


0.03 


0.01 


-0.03 


0.02 


[Ca I/Fe] 


0.01 


-0.02 


0.00 


0.01 


[Sc II/Fe] 


-0.02 


-0.02 


0.00 


-0.05 


[Ti I/Fe] 


0.01 


-0.03 


-0.03 


-0.03 


[Ti II/Fe] 


-0.02 


-0.01 


O.OI 


-0.03 


[Cr I/Fe] 


0.01 


-0.02 


-0.03 


-0.02 


[Mn I/Fe] 


0.01 


-0.02 


-0.04 


-0.03 


[Fe I/Fe] 


0.02 


0.01 


-0.03 


0.01 


[Fe II/Fe] 


-0.03 


-0.02 


0.03 


-0.02 


[Co I/Fe] 


0.01 


-0.03 


-0.04 


-0.03 


[Ni I/Fe] 


0.01 


-0.01 


-0.04 


-0.03 


[Sr II/Fe] 


-0.02 


0.01 


-0.01 


-0.04 


[Ba II/Fe] 


-0.02 


0.00 


-0.01 


-0.04 



dex). Even if we use exactly the same lines in the giants as in 
the dwarfs, we observe the same effect; thus, an error in log 
gf values cannot explain the difference. On the other hand, to 
reduce the derived [Ti/Fe] by 0.2 dex (keeping the same tem- 
perature) would require changing log g in the turnoff stars by 
about 1 dex, which is quite incompatible with the ionisation 
equilibrium of the iron lines. 



7. The light odd-Z metals: Na, Al, K, and Sc 

7. 1 . Sodium and Aluminium 

In both dwarf and giant EMP stars, Na and Al abundances 
can only be derived from the resonance lines, which are very 
sensitive to NLTE effects (Cayrel et al. '2004'). The Na and 
Al abundances in our two stellar samples have been derived 
using the NLTE line formation theory by Andrievsky et al. 
( |2007j and Andrievsky et al. (2008) for Na and Al, respec- 
tively. When NLTE effects are taken into account, the [Na/Fe] 
and [ Al/Fe] abundance ratios are found to be constant and equal 
in the dwarfs and giants in the interval -3.7 < [Fe/H] < -2.5 
([Na/Fe] = -0.2 and [Al/Fe] = -0.1). This can be appreciated 
visually by looking at figure 7 of Andrievski et al. (2007) and 
figure 3 of Andrievski et al. ( I2008I I. 




Si 



mi-. . 




[Fe/H] 

Fig. 4. [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe] in our program 
stars. Symbols as in Fig.|2] 



7.2. Potassium and Scandium 

The K lines are very weak in our EMP turnoff stars, and [K/Fe] 
could not be determined. 

The Sc abundance in the dwarf stars has been measured 
from the Sc II line at 424.6 nm. In giants, 7 Sc lines could be 



used, and the scatter in the abundances from individual lines is 
very small (below 0.1 dex). There is a systematic difference of 
about 0.2 dex between the Sc abundances in the giants and the 
dwarfs (Fig.|5]l. 
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8. Iron-peak elements 

8.1. Chromium, Cobalt, and Nickel 

Fig. |6] shows the [Cr/Fe], [Co/Fe], and [Ni/Fe] ratios for our 
dwarf and giant samples. There is rather good agreement for 
Ni, but [Cr/Fe] and [Co/Fe] are about 0.2 dex higher in the 
dwarfs than in the giants. Recently Lai et al. (2008) have 
measured the chromium abundance in a sample of giants and 
turnoff stars in the same range of metallicity. The same shift 
appears between their giants and turnoff stars (Fig.|7]i. 

Lai et al. also found an offset between the abundances de- 
rived from Cr I and Cr II. Cr II can only be measured in giants, 
and only a single Cr II line (A - 455.865nm) appears at the 
edge of our blue spectra, but the same offset as observed by Lai 
et al. is clearly visible in our data (Figure[8ll. 

The discrepancy between Cr I and Cr II, and between gi- 
ants and turnoff stars, may point to non-LTE effects. The main 
Cr I Unes are resonance lines. Unfortunately no precise struc- 
ture model for the Cr atom exists, so it is not possible to explore 
this hypothesis at present. If significant NLTE effects were con- 
firmed, the most reliable abundances should be those from the 
Cr II line, suggesting that [Cr/Fe] ~ -i-O.l at low metallicity. 

Nissen & Schuster ( |1997| l found a close correlation be- 
tween the abundances of Na and Ni in the interval -0.7 < 
[Fe/H] < -L3. To explain this correlation, it has been sug- 
gested that the production of ^^Ni during an SN II event de- 
pends on the neutron excess, which itself depends mainly on 
the amount of ^^Na produced during hydrostatic carbon burn- 
ing. However, this correlation is not observed in our sample 
(Fig.|9]l. In fact, [Ni/Fe] and [Na/Fe] have the same mean value 
in turnoff stars as in unmixed giants ([Ni/Fe] = 0.0 and [Na/Fe] 
= -0.2). [Na/Fe] is larger in several of the mixed giants, but 
this is due to mixing with the H-burning shell in layers that are 
sufficiently deep to bring products of the Ne-Na cycle to the 
surface (see Andrievsky et al. i2007i) . 

8.2. Manganese 

The Mn abundances have been derived by fitting synthetic 
spectra to the observations, taking into account the hyperfine 
structure of the lines. We noted in Paper V that, in the gi- 
ant stars, Mn abundances determined from the resonance lines 
were lower than those from the lines of higher excitation poten- 
tial by about 0.4 dex. At this stage, we prefer the abundances 
from the high-excitation lines, because the resonance lines are 
more susceptible to non-LTE effects. However, in the five most 
metal-poor giants only the resonance triplet is detected, so for 
these stars the Mn abundance was determined from the triplet 
and corrected by the adopted 0.4 dex offset. 

For most of the turnoff stars analysed here, again only the 
resonance triplet can be detected. In Fig. [TOh the Mn abun- 
dances from these lines have been systematically increased by 
0.4 dex, while in Fig. [TOb [Mn/Fe] is derived from the reso- 
nance triplet profiles in all the stars and plotted without any 
correction. 

In both cases we find a systematic abundance difference of 
about 0.2 dex between the giants and the dwarfs. 




-3.5 -3 
[Fe/H] 

Fig. 5. [Sc/Fe] in the programme stars. Symbols as in Fig.|2] 



8.3. Zinc 

Zinc cannot be unambiguously assigned to the iron-peak cate- 
gory, since it may be formed by a-rich freeze-out and neutron 
capture as well as by burning in nuclear statistical equilibrium. 
In our sample, the only usable line is the strongest Zn i line of 
Mult. 2, at 48 1 nm. The line is very weak in all our stars, and 
we only consider it reliably detected and provide a measure- 
ment when the equivalent width is larger than 0.35 pm. Thus, 
Table [3] gives only six measurements and eleven upper limits; 
for two stars, the spectrum was affected by a defect, and it is 
not even possible to provide an upper limit. 

Figure [TT] shows [Zn/Fe] versus [Fe/H]; upper limits are 
shown as downward arrows and the giant stars from Cayrel 
et al. (2004) as open circles. The upper limits are consistent 
with the trend defined by the giant stars, but the six actual 
measurements appear to define a similar trend, shifted up- 
wards by about 0.4 dex. This could be another example of the 
dwarf/giant discrepancy found for some other elements. 

Since the majority of our [Zn/Fe] data are upper limits, we 
use survival statistics to analyse them. The giant stars with 
[Fe/H]> -3.0 show a constant level of [Zn/Fe]= +0.199 + 
0.080. We selected the dwarf stars in the same metallicity range 
and used asurv Rev 1 . 2 Q (Lavalley et al. 1992) to compute 
the Kaplan-Meier statistics, as described in Feigelson & Nelson 
(119851) . The mean is +0.491 ± 0.055; since the lowest point is 
an upper limit, it has been changed to a detection to compute 
the Kaplan-Meier statistics, which implies that this mean value 
is biased. The comparison of the two mean values, for giants 
and dwarfs suggests that they are only marginally consistent: 
The 75th percentile of the [Zn/Fe] values for dwarfs (+0.223) 
corresponds to the mean value for giants. Changing the upper 
limits to 2cr or 3cr would push the mean value for dwarfs even 
higher, thus making the values of dwarfs and giants even more 
inconsistent. 

We used the generalized version of Kendall's t (Brown et 
al. 119741 ). as described in Isobe et al. ( 1986), to check if there 
is support for a correlation between [Fe/H] and [Zn/Fe] for the 
dwarf stars. The sample is composed of 6 detections and 1 1 
upper limits. The probability of correlation is 91.3%, so there 
is a hint of a correlation, but no conclusive evidence. 



http://astrostatistics.psu.edu/statcodes/asurv 



8 



P. Bonifacio et al.: First Stars XII. Detailed abundances in BMP dwarfs 



For Zn it appears unlikely that the giant/dwarf discrepancy 
is due to NLTE effects: Takeda et al. ( [2005] l computed NLTE 
corrections for the Zn i line at 48 1 nm; the corrections are small 
and negative for metal-poor giants, positive for metal-poor TO 
stars. Thus, if we applied these corrections to our sample, the 
discrepancy would increase from 0.2 to ~ 0.4 dex. 

It is surprising that several stars display upper limits which 
are lower than the [Zn/Fe] ratios found in other stars of simi- 
lar metallicity, suggesting a real cosmic scatter in the Zn abun- 
dance. It is interesting to note that while for giant stars Lai et al. 
(I2OO8I 1 are in good agreement with our determinations, the two 
dwarf stars for which they have Zn measurements appear to be 
in line with the measurements of giants. This may give further 
support to the idea of a cosmic scatter of Zn abundances, or to 
the existence of a Zn-rich population. 

However, it should be kept in mind that the available Zn 
lines are all very weak (detections are about 0.4 pm, upper lim- 
its 0.1-0.2 pm), and the data should not be overinterpreted. We 
have, perhaps somewhat naively, placed the upper limit at the 
measured value for all stars below our chosen threshold. Had 
we decided to put the upper limit at 3cr above the measured 
EW, all the upper limits would move up among the measure- 
ments or beyond, and there would be no hint of any scatter in 
the Zn abundance. From the point of view of survival statis- 
tics, the fact that the standard deviation from the mean is small, 
compared to observational errors, does not support the pres- 
ence of a real dispersion. The question of a scatter in Zn abun- 
dance in EMP dwarf stars clearly needs further study, if possi- 
ble based on different lines. 



9. Neutron-capture elements 

Very few neutron-capture elements can be measured in turnofF 
stars, because their lines are generally very weak. We could, 
however, measure Sr abundances from the blue resonance line 
of Sr II, and sometimes also Ba abundances from the Ba II line 
at 455.4 nm. The Ba line is generally weak (about 0.5 pm) and 
located at the very end of the blue spectrum, where the noise 
is higher As Fig. [12] shows, we find good agreement between 
dwarfs and giants, although the star-to-star scatter is very large, 
as has already been observed for the giant stars. In Fig. [TSlwe 
show the [Sr/Ba] ratio as a function of [Ba/H]. As already no- 
ticed in Paper VIII (Fig. 15), the scatter in this plane is greatly 
reduced. The dwarf stars appear to behave exactly in the same 
way as giants. 

We have recently studied the Ba abundance in dwarfs and 
giants taking into account non-LTE effects (Andrievsky et al. 
I2009I I. but after correcting for NLTE the general behaviour of 
this element remains the same. 



10. Comparison with other investigations. 

Several other groups have now published detailed analyses of 
EMP stars similar to our own, and it is interesting to compare 
their results to ours. We focus on the results of the OZ project 
(Cohen et al. 2004, 2008) and Lai et al. (2008). The details of 
the comparison are provided in appendices |A] |B] and |C] The 
final conclusion of this comparison is that there is excellent 
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Fig. 6. [Cr/Fe], [Mn/Fe], [Co/Fe], and [Ni/Fe] in the program 
stars. Symbols as in Fig.|2] 



agreement between the three groups, and the small differences 
can be understood in terms of differences in the adopted at- 
mospheric parameters, model atmospheres or line selection. 
The MARCS model atmospheres used by us agree with the 
ATLAS non-overshooting models adopted by Lai et al. (I2008I I. 
and both yield abundances which are about 0. 1 dex lower than 
the ATLAS overshooting models adopted by the OZ project. 



11. Discussion 

For most elements, the overall abundance trends defined by 
dwarfs and giants show good agreement. For example, [Ti/Fe] 
is constant at low metallicity in both giants and dwarfs, and 
[Mn/Fe] decreases with metallicity in both giants and dwarfs. 

However, some elements show systematic shifts in [X/Fe] 
between turnofF stars and giants of the same metallicity. 
Generally, [X/Fe](j„arfs - [X/Fe]giants ~ +0.2 dex, except for 
Mg and Si, which show a negative shift. Also, [Cr/Fe] appears 
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Fig. 7. [Cr/Fe] in our stars (circles) compared to those of Lai et 
al. (2008, triangles). Filled symbols: turnoff stars; open sym- 
bols: giants. 





Cr 











-4 -3.5 -3 -2.5 -2 

[Fe/H] 

Fig. 8. [Cr/Fe] in our stars from Cr I lines (open circles: gi- 
ants, filled circles: turnoff stars), and from the single Cr II line 
(crosses; giants only). The large offset cannot be explained by 
measurement errors. 



to be flat in the dwarfs, but displays a significant slope for the 
giants. It is difficult to explain these shifts by systematic errors 
in the models (error in temperature or in gravity) because the 
effects on the abundance of all the elements are very similar 
(see Table |4ll, so the ratios [X/Fe] are little affected. 

These differences are rather puzzling because, except for 
C, N, and possibly Na, the chemical composition of the gi- 
ant stars should be unaltered since the star formed, so one 
would expect that the abundances in giants should match those 
in dwarfs at any given metallicity. The discrepancy we find is 
most likely due to shortcomings in our analysis, but we do not 
know whether we should trust the derived results for giants or 
dwarfs (or perhaps neither!). 

In the following we discuss the two main simplifications of 
standard model atmospheres, the neglect of effects of granu- 
lation ("3D effects" for short) and deviations from local ther- 
modynamic equilibrium (NLTE), as possible causes of the ob- 
served discrepancies. 

11.1. Granulation (3D) effects 

It is well known that hydrodynamical simulations ("3D mod- 
els") predict much cooler temperatures in the outer layers of 



metal-poor stars than ID models (Asplund et al. 1 1 9991 Collet 
et al. 12007' Caffau & Ludwig 120071 Gonzalez Hernandez et 
al. |2008 Paper XI). The effect is more pronounced for dwarfs 
than for giants. The species most affected by this difference 
are clearly those which predominantly reside in such cool lay- 
ers, most notably the diatomic molecules such as CH and NH. 
Since one of the most striking differences between dwarfs and 
giants is in fact the C abundance, which we derive from CH 
lines, we decided to investigate the effects of granulation in 
more detail. 

To accomplish this, we used the two CO^BOLD (Freytag 
et al. 2002, Wedemeyer et al. 2004) 3D models described in 
Paper XI (Tes/log g/[Fe/H]: 6550/4.50/-3.0 and 5920/4.50/- 
3.0). Unfortunately, we do not yet have any fully relaxed mod- 
els for giant stars, so we decided to use a representative snap- 
shot of a 3D simulation of a giant close to relaxation (Teft/log 
g/[Fe/H]:4880/2.00/-3.0). 

Table |6]Usts the mean 3D corrections as defined by Caffau 
& Ludwig ( 120071 ) for the three models described above. The 
sense of the correction is always 3D-1D. Approximating the 
3D correction for the G-band as the average for just 4 lines 
is admittedly somewhat crude, but should provide a reliable 
order-of-magnitude estimate for the effect. 

For the C abundance, the effect is quite prominent for 
dwarfs. The magnitude of the correction is such that, if applied, 
the discrepancy in [C/Fe] between dwarfs and giants would be 
somewhat reduced (from 0.27 dex to 0.13 dex), but with the 
opposite sign, the dwarfs now showing a slightly lower C abun- 
dance. Given the crudeness of our 3D computations, we cannot 
claim with certainty that 3D effects will explain the discrep- 




Fig. 9. [Ni/Fe] vs. [Na/Fe] in dwarfs and giants; symbols as in 
Fig.|2] The dashed line shows the correlation found by Nissen 
& Schuster 09971) for -0.7 < [Fe/H] < -1.3 and corresponds 
to the expected production ratios of Na and Ni in Type II super- 
novae. We do not observe this correlation in our sample. The 
few large [Na/Fe] values ([Na/Fe] > -i-O.l) refer to some of the 
more extreme "mixed" giants discussed in Paper IX. 
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ancy. To the extent that our order-of-magnitude estimates are 
reliable, it is possible that more accurate computations with a 
larger set of parameters, encompassing the full range of our 
dwarf and giant stars, would yield [C/Fe]~ 0.2 for both dwarfs 
and giants. 

For the giant model, our computed correction for C is a fac- 
tor of two smaller than the results of Collet et al. (I2007I I. Also, 
for Fe, our corrections are considerably smaller than found by 
Collet et al., especially for the resonance line. 

The issue clearly requires further investigation, which we 
will undertake when we have several fully relaxed 3D mod- 
els of giants. A detailed discussion is therefore premature. At 
present it is unclear whether the different results we find are 
due to some fundamental difference between the 3D codes: 
CO^BOLD (Freytag et al. 12002 ^, Wedemeyer et al. l2004l l in 
our case, and the Stein & Nordlund (1998 ) code in the case 
of Collet et al. (2007), or simply to the c hoice of a MARCS 
model as the ID reference by Collet et al. (! 2007| l. 

We performed spectrum synthesis computations, using 
Linfor3EQ to estimate 3D corrections for a few selected lines. 



Table 5. Lines used to test the granulation effects. 



http://www.aip.de/ mst/Linfor3D/linfor_3D_manual.pdf 



Species 


A 


X 




nm 


eV 


CH 


430.0317 


0.00 


CH 


430.0587 


0.36 


CH 


430.1072 


1.44 


CH 


430.1135 


0.31 


Si I 


390.5523 


1.91 


Si I 


410.2916 


1.91 


Sen 


424.6822 


0.31 


Tin 


376.1323 


0.57 


Tin 


391.3468 


1.12 


Cri 


425.4332 


0.00 


Cri 


520.8419 


0.94 


Mm 


403.0753 


0.00 


Mm 


404.1355 


2.11 


Fe I 


382.4444 


0.00 


Fe I 


400.5242 


1.56 


Fe I 


418.7795 


2.42 


Fe I 


422.7427 


3.33 


Co I 


384.5461 


0.92 


Zni 


481.0528 


4.08 





-4 -3.5 -3 -2.5 -2 

[Fe/H] 

Fig. 10. [Mn/Fe] ratios for our dwarf and giant stars. In panel 
a), [Mn/Fe] in most of the giant stars is determined only from 
lines with excitation potential > 2eV (open squares). In the 
turnoff stars (black dots) and the five most metal-poor gi- 
ants (open circles), only the resonance triplet of Mn is usable; 
[Mn/Fe] is then derived from these lines, corrected by -1-0.4 
dex (see text). In panel b), [Mn/Fe] is determined from the res- 
onance lines for all stars (giants and dwarfs), without any cor- 
rection. In both cases, the mean [Mn/Fe] ratio is offset by ~0.2 
dex between the giant and dwarf stars. 



listed in Table |5] This is not meant to substitute for a full 3D 
investigation of the sample, but should provide an indication 
whether the differences found between giants and dwarfs might 
vanish if suitable 3D models were used. 

For Si, Co, and Zn the predicted corrections are the same 
for giants and dwarfs, so the discrepancy for these elements 
should not be due to granulation effects. For Sc and Ti, how- 
ever, the differences go in the direction of increasing the dis- 
crepancy between dwarfs and giants. 

For Mn and Cr, the difference in correction between dwarfs 
and giants is such as to exactly cancel the discrepancies. We 
caution, however, that the corrections listed in Table |6] are the 
average of those for the resonance and high excitation lines. 
The difference in correction between the two lines is smaller 
for the giant model (0. 1 dex for Cr, 0.3 dex for Mn) than for the 
dwarf models (0.4-0.5 dex for Cr, 0.7 dex for Mn). This differ- 
ence is still somewhat problematic, however, in the sense that 
while the ID analysis achieved a good excitation equilibrium 
for Cr, an analysis based on the 3D atmospheres does not. This 
suggests that the temperature scale appropriate for 3D models 
may in fact be different from those adopted in this paper and 
by Cayrel et al. (120041 ). As mentioned above, a ID analysis im- 
plies a Mn abundance about 0.4 dex lower for the resonance 
lines than for the high-excitation lines, and the 3D corrections 
for Mn increase this difference, up to 1.1 dex. 

It is unlikely that the use of 3D models will bring the abun- 
dances in giants and dwarfs into agreement for all elements, 
although it may be possible for a few (most likely C, Cr, and 
Mn). However, a full re-analysis based on 3D models, includ- 
ing a redetermination of the atmospheric parameters, is needed 
before reaching a firm conclusion on this point. For the time be- 
ing, since the predicted 3D corrections are always smaller for 
our giant model than for the dwarf models, we consider the ID 
abundances for giants to be more reliable than for the dwarfs. 
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Table 6. Mean 3D corrections for selected elements. 
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Fig. 11. [Zn/Fe] ratios in dwarf stars (this paper; filled circles) 
and in giants (Paper V; open circles). 



1 1.2. Deviations from local thermodynamic equilibrium. 



The analysis in this paper and in Cayrel et al. (I2004I I is based 
on the assumption of local thermodynamic equilibrium (LTE), 
both in the computation of the model atmospheres and in 
the line transfer computations. For Na and Al, results based 
on NLTE line transfer computations have been presented in 
Andrievsky et al. ( |2007l |2D08j . For both elements, the LTE 
computations implied a discrepancy between dwarfs and gi- 
ants, while the NLTE computations provided consistent abun- 
dances between the two sets of stars. In the case of Na, the 
NLTE corrections are not very different for dwarfs or giant 
models for lines of a given equivalent width, but the cor- 
rection depends strongly on the equivalent width. The giant 
stars, which are cooler, have larger equivalent widths and larger 
NLTE corrections. In this case the LTE abundances of dwarfs 
are to be considered more reliable than those of giants. 

The result cannot be generalized, however, so detailed 
NLTE computations should be carried out for all the elements 
for which we find a discrepancy between dwarfs and giants. 
Also, from the point of view of departures from NLTE, one 
cannot a priori assume that the departures are larger for the 
stronger lines (i.e. for giants), although this is often the case. 

Accordingly, except for the two elements Na and Al for 
which we already have NLTE computations, we cannot at 
present say whether accounting for NLTE effects could remove 
the discrepancy between dwarfs and giants. Computations of 
the NLTE abundance of Mg are under way. 

11.3. Could the dwarf/giant discrepancy be real ? 

For C, the difference in [C/Fe] between dwarfs and giants might 
represent the effect of the first dredge-up, which could be re- 
sponsible for a decrease of the C abundance due to a first mix- 
ing with the H-burning layer, where C is transformed into N. 
For the other elements we see no possible nucleosynthetic ori- 
gin for the dwarf/giant discrepancy. 



Another possibility is that the abundances in EMP turnoff 
stars are seriously affected by diffusion (see e.g. Korn et al. 
120071 and Lind et al. 2008 ). From Table 2 of Lind et al. (|2008] | 
one can deduce the following variations in abundance ratios 
between TO stars and RGB stars in the globular cluster NGC 
6397: A[Mg/Fe] = -0.04 + 0.17, A[Ca/Fe] = +0.06 + 0.13, 
A[Ti/Fe] = +0.16 + 0.12. So only for [Ti/Fe] is a variation 
marginally detected, which happens to be of the same order of 
magnitude and sign as the giant/dwarf discrepancy observed by 
us. 

Although a role of diffusion cannot be ruled out, the evi- 
dence in favour is, at best, very weak. Confirmation of the re- 
sults of Korn et al. (2007) and Lind et al. (2008) by an inde- 
pendent analysis would be useful, especially in view of the fact 
that previous investigations of the same cluster (Castilho et al. 
120001 Gratton et al. 2001 ) gave different results. As we pointed 
out in Paper Vll, the adoption of a higher effective tempera- 
ture for the turn-off stars of this cluster, as done by Bonifacio 
et al. (I2OO2I 1. would largely cancel the abundance differences 
between TO and RGB. Even if the results for NGC 6397 were 
confirmed, it is not obvious that they would apply to the field 
stars analysed in the present paper Unlike the stars in a glob- 
ular cluster, these stars are not necessarily strictly coeval, and 
their metallicities range from ~0.7 to ~1.7 dex below that of 
NGC 6397. 

1 1.4. Do the giant and dwarf samples belong to the same 
population? 

It could be argued that the observed giant and turnoff sam- 
ples might belong to different populations, since the giants 
would, on average, be more distant than the turnoff stars. To 
test this, we have compared the radial velocities of the two 
samples (it would have been preferable to compare the space 
velocities, but the distances and proper motions of the giants 
are generally very uncertain). Barycentric radial velocities for 
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the turnoff stars are given in Bonifacio et al. (120071 ). For the 
giants they are given in Table [Tj they are based on the yellow 
spectra centered at 573nm with laboratory and measured wave- 
lengths of numerous Fe i lines (Nave et al. 119941 ). The wave- 
lengths for the telluric lines for the zero points have been taken 
from Jacquinet-Husson et al. (I2005I I. Velocity errors should be 
below 0.3 km s \ more than adequate for the present purpose 
(see also Hill et al. 2002 ). 

Since all the program stars (except for a few of the giants) 
have been selected from the HK survey (Beers et al.l l985|[1992| 
and Beers 119991 ). which is kinematically unbiased, their ra- 
dial velocities should be an unbiased estimate of the kinematic 
properties of the population. Thus, if the stars were indeed 
drawn from different populations, we would expect their radial- 
velocity distributions to differ The mean radial velocities and 
standard deviations are -12 and 141 kms^'for the giants, - 
32 kms"' and 159 km s"' for the turnoff stars, respectively. 
A Kolmogorov-Smirnov test shows only a 10-15% probability 
that the two samples have not been drawn from the same parent 
population. Thus, the radial-velocity data support the assump- 
tion that the dwarfs and giants belong to the same population. 

12. Conclusions 

We have determined abundances of C, Mg, Si, Ca, Sc, Ti, Cr, 
Mn, Co, Ni, Zn, Sr and Ba for a sample of 18 BMP turnoff stars, 
which complements the sample of giants discussed by Cayrel 
et al. (12004'). For the subgiant BS 16076-006 it was possible 
also to determine the N abundance. 

For Ca, Ni, Sr, and Ba we find excellent consistency be- 
tween the abundances in dwarfs and giants at any given metal- 
licity. For the other elements we find abundances for the dwarfs 
which are about 0.2 dex larger than for giants, except for Mg 
and Si, for which the abundance in dwarfs is about 0.2 dex 
lower than in the giants, and Zn, for which the abundances in 
dwarfs are about 0.4 dex higher than in the giants. 

The only element for which such a discrepancy could have 
an astrophysical explanation is C. In fact, if the first dredge-up 
were capable of bringing into the atmosphere material which 
had undergone CN processing, one would expect to find lower 
C abundance in giants than in dwarfs. Such an effect is not 
predicted by standard models of stellar evolution and would 
require some extra-mixing mechanism. For all the other ele- 
ments which display a discrepancy between dwarfs and giants 

T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-^ 

O 1 

o 6>(?§, ° 'I. i§> ■ 
O-qO '..Q-A. o^....o _ 

o o o ° n : 

o o Q . 

-6 -5 -4 -3 -2 -1 

[Ba/H] 

Fig. 12. [Sr/Ba] as a function of [Ba/H] in our dwarf (dots) and 
giant stars (open circles; data from Paper VIII). 
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Fig. 13. [Sr/Fe] and [Ba/Fe] ratios in our dwarf and giant stars; 
symbols as in Fig. 121 Note the large scatter: The vertical scale 
is not the same as in the other figures. 



we are unable to find any plausible astrophysical explanation. 
We conclude that the discrepancies arise from shortcomings in 
our analysis, probably also for C, but certainly for all other el- 
ements for which discrepancies are found. 

We have made an approximate assessment of the efi'ects of 
granulation and conclude that they are unlikely to explain the 
discrepancies, except perhaps for C, Mn and Cr. In any case, the 
3D corrections appear to be smaller for giants than for dwarfs, 
which suggests that the ID abundances of giants are preferable 
as reference data for studies of the chemical evolution of the 
Galaxy. 

The other obvious shortcoming in our analysis is the as- 
sumption of local thermodynamic equilibrium. Detailed NLTE 
line transfer computations exist for Na and Al (Andrievsky et 
al. 120071 12008J . and for these two elements they in fact re- 
move the dwarf/giant discrepancy implied by the LTE analy- 
sis. Computations for Mg are in progress, and it seems that 
the agreement between giants and dwarfs is at least improved. 
This result cannot be generalized to other elements, and it is 
not clear whether NLTE computations might remove any of 



the other discrepancies. Clearly, NLTE computations for other 
key elements are urgently needed. 

For readers who wish to use our data for comparison with 
Galactic evolution models we suggest that, for elements for 
which a dwarf/giant discrepancy exists, the abundances in gi- 
ants are to be preferred. We plan to publish an updated table 
of all the abundances in the First Stars programme in a final 
paper of the series. For the time being we direct the reader who 
wants the most updated abundances of the First Stars giants, 
to the following papers: for Li, C, N and O, Spite et al. (120051 
First S tars VI; ^006, First Stars IX); for Na, Andrievsky et al. 
( I2007D ; for Mg to the NLTE abundances in Fig.|4]of the present 
paper (to be published in full soon); for Al, Andrievsky et al. 
( I2UO8I 1; for K, Ca, Sc, Ti, Mn, Fe, Co, Ni, Zn, Cayrel et al. 
( 120041 ); for Cr, the Cm abundances given in Fig. [8] should be 
preferred; for Ba Andrievsky et al ( |2009l l; for all the other ele- 
ments heavier than Zn Frangois et al. ( 120071 First Stars VIII). 

The reasons for recommending the use of abundances in gi- 
ants are threefold: 1) granulation effects are smaller for giants 
than for dwarfs; 2) giants have lower effective temperatures and 
stronger lines, so from the observational point of view their 
abundances are better determined; and 3) the atmospheres of 
giants stars are well mixed by convection and should be im- 
mune to chemical anomalies driven by diffusion. One should, 
however, bear in mind that future NLTE analyses of our data 
could imply substantial revision of the abundances in both gi- 
ants and dwarfs. 
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Fig.A.l. Temperature structure for three models with Teff= 6365, logg4.4 and [M/H]= -3.0. The sohd line is our MARCS 
models, the dashed line is an atlas overshooting model, the dashed-dotted line is an ATLAS non-overshooting model. 



Appendix A: Comparison with the OZ project. 

The OZ project (Cohen et al. 120041 l2008l l has produced a data set similar to that of the "First Stars" project, it is therefore of some 
interest to verify how these data sets compare. Cohen et al. (2004) analysed a set of dwarf stars which is directly comparable to 
those analysed in the present paper. The spectra were acquired with the HIRES spectrograph at the Keck I Telescope, at a resolu- 
tion only slightly lower than our UVES-VLT data (34000 rather than 45 000), and the S/N ratios are comparable. The equivalent 
widths were measured using an automatic code which fits gaussians, therefore the general philosphy of EW measurement does 
not differ from ours. In fact Cohen et al. (I2008I I observed two giant stars measured by Cayrel et al. (I2004I I. and the equivalent 
widths compare very well (see figures 13 and 14 of Cohen et al. 120081 and related text). The two projects differ in the method 
used to fix the atmospheric parameters: we use the wings of Ha for dwarf stars, while the OZ project relies on photometry to 
derive Tefj. For surface gravity we use the iron ionisation equilibrium, while the OZ project relies on theoretical isochrones. 

We have also investigated the gf values used by the two projects, and they are very similar; the use of the one or the other set 
would not imply differences in the derived abundances smaller or equal to 0.02 dex. 

Thus, part of the differences will depend on the different adopted atmospheric parameters. There is no dwarf star in common 
between the two groups; thus it is not straightforward to compare the results of the two projects. 

For the analysis the two projects use different model atmospheres and different line formation codes. We use MARCS model 
atmospheres and turbospectrum, while the OZ project uses ATLAS models interpolated in the grid of Kurucz (1993), with the 
overshooting option switched on, and the MOOG code (Sneden fl 9731 [1974, 2007). As we shall show below, the different choice 
of line formation is relatively unimportant, implying differences in the abundances of a few hundredths of dex; on the other hand 
the choice of ATLAS overshooting models implies abundances which are higher by about 0.1 dex for all the models. Such a 
behaviour was already noticed by Molaro et al. ( 119951 ) for Li, but we show here that it is indeed true for all species. 

In Table lA.ll we list the abundances for the star HE 0508-1555 derived by using the equivalent widths of Cohen et al. ( |2004| l 
and their atmospheric parameters (T^s- 6365, log g - 4.4 and a microturbulent velocity of 1.6 kms"') with three different 
models: a MARCS model interpolated in our grid, an ATLAS model computed without overshooting and an ATLAS model 
computed with overshooting. For all the models we assumed [M/H]=-3.0. Our ATLAS models are somewhat different from 
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Fig. A.2. Temperature structure of the deepest layres of three models with Teff= 6365, \ogg= 4.4 and [M/H]= -3.0. The solid line 
is our MARCS models, the dashed-dotted line is an ATLAS non-overshooting model with q-mlt^ 1-25 (also shown in Fig. IA.lt , 
the dotted line is an ATLAS non-overshooting model with cmlt^ 1 00. 



those of the Kurucz (1993) used by the OZ project. In the first place we use the "NEW" opacity distribution functions (Castelh 
& Kurucz |2003l l computed with 1 km s ' microturbulence. In the second place we use the Linux version of ATLAS (Sbordone 
et al. l2004l l. In all cases the line formation code used was turbospectrum. In the last two columns of Table lA.ll we provide the 
abundances of Cohen et al. (2004 ), for the reader's convenience. 

Inspection of Table [ATI immediately suggests that both the difference in ATLAS versions and the different line formation 
codes used are immaterial, since the abundances we find for almost all elements are within 0.04 dex of those of Cohen et al. 
( I2004I I. The two exceptions are Al and Si. For Al there is a good reason for the discrepancy: Both Al i lines used are affected by 
the neighbouring Balmer lines. In our analysis we used spectrum synthesis to derive the abundances. Instead, MOOG can take 
into account the absorption due to the Balmer lines, either using the opacit switch to introduce a fudge factor on the continuum 
opacity or using the strong keyword to read strong lines to be considered. 

For Si the difference between our result with the ATLAS overshooting model and the published value of Cohen et al. (|2004| l 
is 0.08 dex. This abundance is based on a single line of about 10 pm of EW, therefore clearly saturated. The precise value of the 
damping constants used for this line and the way the different codes use them may have an impact. 

Another inference which can be drawn from Table lA.fl is that MARCS models and ATLAS non-overshooting models provide 
results which are quite similar. This is not the case for the ATLAS overshooting models, which imply abundances which are higher 
by about 0.1 dex for all elements. The reason for this behaviour may be understood by looking at the temperature structure of 
the different models. In Fig. lA.ll we compare the temperature structures of our MARCS model (solid line), the ATLAS non- 
overshooting model (dashed-dotted line) and the ATLAS overshooting model (dashed line). The temperature structure of the 
ATLAS non-overshooting and of the MARCS model are quite similar In fact, the only difference is for the deepest layers and is 
driven by the different choice made for the mixing length. 

In Fig. IA.2I we show the temperature structure of the deeper layers of our MARCS model (solid line) together with two 
ATLAS non-overshooting models with different values of qtmlt 1-25 (dashed-dotted line) and 1.00 (dotted line). The ATLAS 
model with Q'mlt= 1 -00 is closer to the MARCS model, up to log T500 ~ 0.7, but then becomes hotter than the MARCS model. In 
general it is impossible to chose a ckmlt such as a MARCS and an ATLAS model have exactly the same structure in depth, due 
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to the different formulation of the mixing-length therory in the two codes. Such differences in the very deepest layers have very 
little influence on a typical abundance analysis. In fact only the lines which form in these very deep layers are affected, i.e. very 
weak lines of 0. 1 pm or smaller, and the wings of Hy6 and higher members of the Balmer series. 

In general we can conclude that MARCS and ATLAS non-overshooting models are very similar, and an abundance analysis 
based on the two models will yield abundances which are consistent within a few hundredths of dex. The situation is dramatically 
different when we consider the ATLAS overshooting models. Such models present a temperature structure which is very different 
from both ATLAS non-overshooting and MARCS models in the region -1 < logrsoo < 1 where the majority of lines used in 
abundance analysis are formed. 

Castelli, Gratton & Kurucz ( |1997al |1997bl l investigated extensively the effects of the approximate overshooting present in 
ATLAS and concluded that the no-overshooting models are capable of reproducing a larger set of observables, thus discouraging 
the use of overshooting models. To these considerations we may add that having investigated the mean temperature structures of 
CO^BOLD 3D hydrodynamical models we never see the typical "bump" in the temperature structure seen in ATLAS overshooting 
models. The real effect of the overshooting is the over-cooling of the outer layers with respect to what is predicted in radiative 
equiUbrium models (Asplund et al. [19991 Collet et al. l2007l Caffau & Ludwig 2007 Gonzalez Hernandez et al. l2008l Paper XI). 
This is a further reason to avoid the use of the ATLAS overshooting models. 

It can be appreciated that the differences due to different models largely cancel out when considering abundance ratios, like 
e.g. [Mg/Fe], rather than abundances relative to hydrogen. For example [Mg/H] is -2.04 for the ATLAS non-overshooting model, 
but -1.90 for the ATLAS overshooting one, however [Mg/Fe] is 0.50 in the first case and 0.51 in the second case. 

A difference in the average [Mg/Fe] is found between us and the OZ project, of the order of 0.2 dex (the OZ project being 
higher), both if we consider only dwarf stars, only giants, or the full samples. Such an offset is roughly compatible with a Icr 
error on each side, but perhaps a little disturbing. Only a 0.01 dex difference is due to the different adopted solar abundances. The 
use of different models and different atmospheric parameters should largely cancel out when considering a ratio such as [Mg/Fe]. 
Largely does not mean totally, however: Table IC.fl shows a 0.06 dex difference in [Mg/Fe] for BS 16467-062, depending on the 
adopted atmospheric parameters. 

Table 10 of Cohen et al. (2004) is also illuminating; it shows how the average [Mg/Fe] changes if one considers the mean 
computed from the abundances derived from a single line of Mg i. Of the five Mg i lines used by Cohen et al. (2004) three tend 
to give systematically higher abundances, while two give systematically lower abundances. The final result depends on the set 
of adopted lines. This issue requires further investigation in the light of the study of deviations from thermodynamic equilibrium 
for the Mg i lines. Our abundance ratios are in agreement with those provided by the OZ project, within the stated errors. 

At the end of this exercise we conclude that our measurements and those of the OZ Team are highly consistent. Differences 
in the published abundances can be traced back to the different atmospheric parameters adopted, the different treatment of 
convection in the adopted model atmospheres (approximate overshooting versus no overshooting), and for some elements to the 
paiticular choice of hnes. 



Appendix B: Details of the comparison witli Lai et al. 2008 

Lai et al. (2008) also analysed a set of stars which is comparable to that of the First Stars project with respect to metallicity. Their 
sample is also extracted from the HK survey and comprises both dwarfs and giants. Their method to determine atmospheric 
parameters is similar to that of the OZ project, photometric temperatures from the V - K colour and gravities derived from 
isochrones. They observe the giant star BS 16467-062, also observed by us (Paper V) and in the OZ project (Cohen et al. 2008) 
and, not surprisingly, derive atmospheric parameters very close to those of Cohen et al. This allows a very tight comparison of 
the analysis by the three groups, which we defer to Sec.lCl 

Lai et al. use the same spectrum synthesis code as we and also use ATLAS 9 non-overshooting models which, as discussed 
in Sec.|A] are very similar to our MARCS models. It is therefore to be expected that the abundance ratios determined by the two 
groups are quite similar In Fig. lB.ll we compare the [Mg/Fe] ratios of the First Stars project with those of Lai et al. The overall 
agreement is satisfactory. 

In Fig. |B.2| we compare the [O/Fe] ratios of the First Stars project (only giants) with those of Lai et al. The figure seems to 
indicate a good agreement; however we believe that this agreement is in fact fortuitous, as our oxygen abundances were based 
on the 630 nm [OI] line, while those of Lai et al. have been derived from one OH line of the UV A~Y. - X^Yl electronic system 
around 318.5 nm (although the precise line used is not specified). These OH lines are known to provide very high [O/Fe] ratios 
when analysed with ID model atmospheres (e.g. Boesgaard et al. 1999, Israelian et al. 2001). Asplund & Garcia Perez (.2001) 
have explained this behaviour as due to overcooling of the outer layers of the stars, caused by the overshooting of the convective 
elements and not properly described by ID model atmospheres. Our own hdyrodynamical computations (Gonzalez Hernandez 
et al. 2008 Paper XI) confirm this interpretation. In view of this fact it is, at first sight, surprising to find that Lai et al. determine 
rather low [O/Fe] ratios from the OH lines. Closer inspection of their analysis reveals, however, that this is mainly driven by their 
adopted g f values for these lines. 
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Fig.B.l. Comparison of the [Mg/Fe] ratios of the First Stars project and those of Lai et al. (2008). Our data is shown as circles, 
while those of Lai et al. as triangles. Open symbols correspond to giant stars, while filled symbols to dwarfs. 



In Fig. IB. 31 we show a portion of the spectrum of CS 31085-024, used by Lai et al., which we downloaded from the Keck 
Observatory Archive^ compared with two synthetic spectra computed using an ATLAS 9 model with the atmospheric parameters 
adopted by Lai et al. and two different OH line lists. In the first case (red line) we adopted the gf values for the OH lines of the 
(0-0) vibrational band of the A^S - X^Yl electronic system computed from the lifetimes of Goldman & Gillis ( I1981I I. which we 
used in Paper IX. In the second case (green spectrum) we used the lines computed by R.L. Kurucz. This second list is far richer, 
since it includes also lines from other vibrational bands and not only the (0-0) band. However, even from this limited portion of 
the spectrum it can be appreciated that the Kurucz gf values are larger than those derived from the Goldman & Gillis ( |1981| l 
lifetimes; use of the latter gf values would lead to considerably larger OH abundances. 

For this reason we believe that the oxygen abundances in the stars of the Lai et al. sample should be reinvestigated using 
a different set of gf values and hydrodynamical model atmospheres. It is likely that the 3D corrections for the giant stars (the 
majority of the Lai et al. sample with oxygen measurements) are smaller than those for dwarf stars (see Paper XI), since the 
overcooling is far less extreme in giants than in dwarfs, it is however unlikely that the effect is negligible. 

We disagree with the statement by Lai et al., who discard the use of 3D models for the analysis of the OH lines since "these 
models seem to overpredict the solar oxygen abundance derived from helioseismology (Delahaye & Pinsonneault 2006)". In the 
first place the oxygen abundance in the Sun is not derived from OH UV lines; in the second place it is now clear that the low 
solar oxygen abundances which have been claimed in the past (Asplund et al. 2004) are not due to the use of 3D hydrodynamical 
models, but to low measured EWs and extreme assumptions on the role of collisions with H atoms in the NLTE computations 
(see Caff'au et al. 2008, for a discussion and a new measurement of the solar oxygen abundance). In our view the use of 3D 
hydrodynamical models is necessary for a reliable analysis of OH lines in metal-poor stars. 

The [Cr/Fe] ratios were compared in Fig. [T] and we see that the picture which emerges is very consistent between the two 
analyses, including the dwarf-giant discrepancy discussed in Sec. 18. II In agreement with us, Lai et al. note that when Cm lines 
are measurable the [Cr ii/Fe] ratio remains close to zero, suggesting that the decrease in [Cr/Fe] with decreasing metallicity, seen 
when Cri lines are used, is probably an artifact due to deviations from LTE. 

Finally in Fig. |B.4| we compare the [Zn/Fe] ratios with those of Lai et al. (2008). The have measured Zn only in two dwarfs, 
slightly more metal-rich than our ones and the Zn abundances for these two are in line with what derived from the giants. We 
note that the gf value adopted by Lai et al. is 0.04 dex lower than adopted by us. 



^ http://www2.keck.hawaii.edu/koa/public/koa.php 
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Fig. B.2. Comparison of the [O/Fe] ratios of the First Stars project and those of Lai et al. (2008). Symbols as in Fig. IB. 11 Our 
oxygen abundances are derived from the 630 nm [OI] line, while those of Lai et al. from one UV OH line of the A^S - X^Tl 
electronic system around 318.5nm. 



Appendix C: Comparison for BS 16467-062 

The giant star BS 16467-062 has been observed independently by all three groups, ourselves (Cayrel et al. 2004, Paper V), 
the OZ project (Cohen et al. 2008) and Lai et al. (2008). The latter two groups have used HIRES @Keck, while we have used 
UVES@VLT. 

In their appendix B, Cohen et al. (2QQE\ make a detailed comparison between the analysis of giant stars analysed by us and 
their own analysis. They conclude that the same star analysed by the two groups will show a difference of 0.3 dex in [Fe/H]. 
This is based on their analysis of the giant star BS 16467-062. We wish to explain how this difference arises. We used the EWs 
of Cohen et al. (120081 ) for this star and their gf values to redetermine the abundances using four models: a MARCS model and 
two ATLAS models (overshooting and non-overshooting) for Teif= 5365 K \ogg= 2.95, which are the parameters of Cohen et 
al. (12008) and the MARCS model with Teif= 5200 K , logg= 2.50, which was used in Cayrel et al. ('2004 Paper V). The results 
are shown in Table IC.ll We omit the results from the ATLAS non-overshooting model, since they are identical to those obtained 
from the MARCS model. This could be expected by looking at Fig. IC.ll in which the temperature structures of the two models 
are compared. 

The differences in the abundances between the MARCS model with the parameters of Paper V and those from an ATLAS 
overshooting model with the higher Tefi and \ogg of Cohen et al. (2008) may indeed be as large as 0.3 dex. However, it is 
important to understand that this difference is due to two distinct factors: on the one hand the change in Tetf and logg, as each 
species displays a slightly different sensitivity to these; on the other hand the use of approximate overshooting in the models of 
Cohen et al. (2008). The two effects are comparable. 

The difference in the abundances obtained from the two different MARCS models allow to estimate the sensitivity of the 
various abundances to the model parameters. The difference between the MARCS and ATLAS overshooting model allow to see 
the effect of the approximate overshooting. We confirm that [Fe/H] for this star is 0.3 dex higher using the parameters of Cohen 
et al. ( I2OO8I 1 and an ATLAS overshooting model, relative to that derived using the parameters of Paper V and a MARCS model 
(or an ATLAS non-overshooting model). However, 0.17 dex of this difference arises from the different choices in Tgff and logg, 
and 0. 13 dex comes about from the use of the approximate overshooting. 

Having understood these differences we may conclude that there is excellent agreement between the two analysis. Note 
that with our MARCS model and atmospheric parameters, but the EWs and gf values of Cohen et al. ( 2008) , [Fe/H] for this 
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Fig.B.3. HIRES-Keck spectram of the dwarf star CS 31085-024. The data is the same used by Lai et al. (2008), downloaded 
from the Keck Observatory Archive (http://www2.keck.hawaii.edu/koa/public/koa.phpV Overimposed on the spectrum are two 
synthetic spectra, computed with SYNTHE, from an ATLAS 9 model with TefF= 5949, \ogg- 4.57 and metallicity -3.0 and 
[0/Fe] = 1.54. The synthetic spectrum plotted in red has been computed using the g f values for the OH lines of the (0-0) vi- 
brational band of the A^Z - X^Yl electronic system computed from the lifetimes of Goldman & Gillis (il981) . Instead the 
one in green has been computed using the line list of the OH A^S - X^Il computed by R.L. Kurucz and distributed through 
( |http://wwwuser.oats.inaf.it/ atmos/tarballs/mol ecules.tar.bz2| ). 



star is -3.80, which compares very well with -3.77 given in Paper V. Note also that when using MARCS models (or ATLAS 
non-overshooting models) our atmospheric parameters achieve a slightly better iron ionisation equilibrium (0.03 dex) than the 
parameters chosen by Cohen et al. (2008 0.06 dex). However since both these difference are much smaller than the line-to-line 
scatter it is impossible to chose which set of parameters is better by just looking at the iron ionisation equilibrium. As noted 
above most of these differences tend to cancel out when considering abundance ratios. 
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Fig. B.4. Comparison of the [Zn/Fe] ratios of the First Stars project and those of Lai et al. (2008). Symbols as in Fig. IB. 11 



Iron is the element for which the largest number of lines is measured and in this respect its abundance is more robust. For 
other elements the difference between the values published in Paper V and an analysis by the OZ Team may also reflect the 
different choice of lines. For instance for BS 16467-062 Cohen et al. (I2008I I measure 4 Mg i lines, while in Paper V we measured 
8 lines, but used only 7 to derive the mean Mg abundance. The Mg lines in BS 16467-062 are all weak; thus, the re-measurement 
of the Mg abundance using line profile fitting (see section ISTl i confirms the abundances provided in Paper V. 

We discarded Mg 1 416.7271 nm because the abundance derived from this line deviates strongly from those derived from the 
other lines. The line is rather weak (0.75 pm as measured in our data or 0.68 pm as measured by Cohen et al. I2008I I. but even 
for these very weak lines the measurements are highly consistent. Thus the mean Mg abundance from our 7 lines is, as given in 
Paper V, 3.97, with a rather small scatter of 0.09 dex. On the other hand the mean Mg abundance from the four lines measured 
by Cohen et al. (l2008l l. including Mg 1 416.7271 nm, and using the atmospheric parameters and model of Paper V, is 4.15 with a 
rather large scatter of 0.33 dex. The mean Mg abundance for these four lines from our measurements is 4.12 with a scatter of 
0.38 dex. Finally if we take the measurements of Cohen et al.(2008^ and discard the Mg i 416.727 1 nm line we obtain 3.99 with 
a scatter of 0. 1 1, highly consistent with our published value in Paper V. 

The three groups (First Stars, OZ project, Lai et al.) have used different atmospheric parameters for this star, and the sensitivity 
of abundances to these is detailed in all the three papers. In order to make a stringent comparison between the results of the three 
groups it is advisable to derive abundances from each set of EWs and gf values for a same model atmosphere and with the 
same spectrum synthesis code. We did so in Table IC.3I in which we used the MARCS model used in Paper V to rederive all 
the abundances. We compare the atomic species in common, excluding Al, for which both ourselves and Lai et al. have used 
spectrum synthesis. 

Inspection of Table |C3] immediatelv reveals that, with very few exceptions, the abundances of the First Stars project rely 
on a larger number of lines than those of the other teams. This is particularly striking for iron, for which we use 130 Fe ii lines 
compared to 55 of Cohen et al. and 52 of Lai et al.; a similar situation is found for Ti, where we use 1 1 Ti i and 23 Ti ii lines while 
Cohen et al. use 2 and 14, respectively and Lai et al. 1 and 9. This probably reflects the fact that the First Stars spectra have a 
larger total wavelength coverage and a more uniform high S/N ratio across the spectra. This is due in part to the fact that UVES, 
as a two-arm spectrograph, covers roughly a 30% larger spectral range in a single exposure than HIRES, and in part to the large 
amount of telescope time invested in the First Stars project. 

Once the Mg i line at 416.7 nm has been removed from the set of Cohen et al. the Mg abundance appears to be in remarkably 
good agreement, in spite of the much larger number of lines used by the First Stars team. 
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Fig. C.l. Temperature structure for three models with Teff= 5365, logg2.95 for BS 16467-062. The solid line is our MARCS 
models, the dashed line is an ATLAS overshooting model, the dashed-dotted line is an ATLAS non-overshooting model. 



That the actual choice of lines does make a difference is obvious if we look at the Ca abundances. There is a difference of 
0.23 dex in the Ca abundance derived in Paper V and that of Cohen et al. (2008). Of the three lines measured by Cohen et al. 
(2008) we have only two. The mean Ca abundance for these two lines is 2.81 with a 0.05 dex deviation, thus the discrepancy is 
reduced to 0.1 dex, totally consistent with the observational eiTors. We have measured all four Ca lines used by Lai et al. (2008), 
and the mean of these four lines is close to the abundance given in Table IC.31 However, Cai 443.5 nm appears to be discrepant 
by 0.39 dex with respect to the mean of the other three lines, which is 2.86, only 0.08 dex higher than the value of Lai et al. and 
fully consistent with observational errors. It is then clear that for the species for which a limited number of lines is available, the 
actual choice of lines can make a difference. 

Another noticeable difference is for Si. All three groups have determined the Si abundance from a single Si line, however 
the other two teams have used the Sii 390.6 nm line, while we have used the 410.3 nm line since the other line is heavily 
contaminated by CH lines in the spectra of giant stars. On the other hand, the EWs for the 390.6 nm line agree well among the 
three investigations (9.18 pm for us, 9.34 pm for Cohen et al. 2008 and 9.06 nm for Lai et al. 2008); thus the Si abundance derived 
from this line agrees well among the three investigations. 

It is reassuring that for iron, for which all three groups have measured a large number of lines, the results are fully consistent. 

The conclusion of these comparisons is that the results of the three teams are consistent, once the different choice of atmo- 
spheric parameters and models has been factored out. Some caution must be exercised for the species which are represented by 
few lines, where the actual choice of lines can make a difference, especially if differential NLTE effects are present. 
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Table 7. Radial velocities of the sample of EMP giant stars. 



n 


Star 


obs. date 


Julian date 


Vr 
km s 


1 


HD2796 


2000-10-13 


51830.1114692 


-60.87 


2 


HD122563 


2000-07-18 


51743.9645395 


-26.39 


3 


HD 186478 


2000-10-12 


51829.9786466 


+31.74 


4 


BD+ 17:3248 


2000-10-19 


51836.9868382 


-146.55 






2001-06-03 


52063.2428859 


-146.45 


5 


BD-18 5550 


2000-10-11 


51828.9894138 


-124.84 






2001-09-05 


52157.0417540 


-124.86 


6 


CD-38 245 


2000-07-19 


51744.3756350 


-1-46.39 






2000-07-19 


51744.3970970 


-1-46.42 


7 


BS16467-062 


2001-06-04 


52064.0495221 


-90.58 






2001-07-04 


52094.9546310 


-90.53 


8 


BS 16477-3 


2001-06-03 


52063.1180688 


-222.61 






2001-06-04 


52064.0946300 


-222.59 


9 


BS 17569-49 


2001-06-02 


52062.3702745 


-213.10 






2001-06-03 


52063.3972078 


-213.13 


10 


CS22169-35 


2000-10-14 


51831.2987792 


-1-14.40 


11 


CS22172-2 


2000-10-18 


51835.2952562 


-1-251.03 


12 


CS22186-25 


2001-10-19 


52201.2882394 


-122.34 






2001-10-21 


52203.3130894 


-122.43 






2001-11-07 


52220.2859610 


-122.35 


13 


CS22189-09 


2000-10-15 


51832.2679605 


-20.22 


14 


CS22873-55 


2001-06-01 


52061.2269695 


+214.54 






2006-10-16* 


54024.9749432 


+214.23 


15 


CS22873-166 


2000-10-19 


51836.0985003 


-17.03 






2006-10-17* 


54025.9801752 


-17.24 


16 


CS22878-101 


2000-07-19 


51744.051286 


-129.12 






2000-07-19 


51744.097628 


-129.24 






2000-08-06 


51762.984927 


-129.32 


17 


CS22885-96 


2000-08-09 


51765.0892940 


-250.28 






2000-08-09 


51765.1072672 


-250.39 






2000-08-11 


51767.0918537 


-250.01 






2000-08-11 


51767.1132987 


-249.97 


18 


CS22891-209 


2000-10-18 


51835.9948297 


+80.29 






2006-10-15* 


54023.9716681 


+79.98 


19 


CS22892-52 


2001-09-08 


52160.1424999 


+ 13.31 


20 


CS22896-154 


2000-10-12 


51829.9948968 


+137.96 






2000-10-14 


51831.9888536 


+137.82 


21 


CS22897-8 


2000-10-12 


51829.0123961 


+266.69 


22 


CS22948-66 


2001-09-05 


52157.1005735 


-170.67 






2006-10-18* 


54026.1255044 


-170.62 


23 


CS22949-37 


(see Depagne et al., 2002) 


-125.64 


24 


CS22952-15 


2000-10-14 


51831.0576878 


-18.24 






2006-10-18* 


54026.1033137 


-18.43 


25 


CS22953-3 


2001-09-08 


52160.3271107 


+208.61 






2001-09-09 


52161.3657380 


+208.44 


26 


CS22956-50 


2001-09-04 


52156.1279571 


-0.09 


27 


CS22966-57 


2001-09-04 


52156.1945760 


+101.22 


28 


CS22968-14 


2000-10-12 


51829.3235079 


+159.16 






2000-10-12 


51829.3662157 


+ 159.09 






2000-10-17 


51834.2971840 


+ 159.16 






2001-11-09 


52222.2098787 


+ 159.02 






2001-10-20 


52202.0858094 


1 AH AH 


30 


CS29495-41 


2001-06-02 


52062.2774334 


+79.61 


31 


CS29502-42 


2000-10-12 


51829.0798142 


-138.05 






2001-09-06 


52158.1633788 


-138.16 






2001-09-09 


52161.1907037 


-138.07 


32 


CS295 16-24 


2001-06-04 


52064.3718149 


-84.27 


33 


CS295 18-51 


2000-10-13 


51830.1969033 


+96.42 


34 


CS30325-94 


2000-08-08 


51764.9714377 


-157.76 






2000-08-08 


51764.9928808 


-157.67 


35 


CS31082-01 


(see Hill et al., 2001) 


+139.07 
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Table A.l. Abundances for HE 0508-1555 for different model atmospheres. 



Ion 


A 


cta 


A 




A 




A (Ta 




MARCS 


ATLAS 


ATLAS 


Cohen 








NOVER 


OVER 


et al. 2004 


Mgi 


5.58 


0.27 


5.54 


0.29 


5.68 


0.27 


5.70 0.28 


All 


3.28 


0.06 


3.24 


0.06 


3.38 


0.06 


3.91 0.09 


Si I 


5.19 




5.14 




5.31 




5.23 


Cai 


4.15 


0.08 


4.11 


0.08 


4.24 


0.08 


4.22 0.11 


Sen 


0.99 


0.03 


0.96 


0.04 


1.04 




1.02 0.03 


Til 


3.18 


0.13 


3.14 


0.13 


3.26 


0.13 


3.23 0.14 


Tin 


3.05 


0.09 


3.02 


0.09 


3.10 


0.09 


3.08 0.10 


Cri 


2.94 


0.03 


2.89 


0.03 


3.01 


0.03 


2.99 0.04 


Mm 


2.31 


0.05 


2.26 


0.05 


2.38 


0.04 


2.35 0.05 


Fei 


5.00 


0.16 


4.96 


0.16 


5.09 


0.15 


5.07 0.17 


Fen 


5.11 


0.16 


5.09 


0.16 


5.16 


0.16 


5.13 0.17 


Co I 


2.88 


0.11 


2.84 


0.11 


2.96 


0.10 


2.92 0.11 


Nil 


3.77 




3.72 




3.86 




3.83 



Table C.l. Abundances for BS 16467-062 different model atmospheres and the EWs of Cohen et al. (2008). 



Ion 


MARCS 


ATLAS 


MARCS 






OVER 






T=5364K 


T=5364 


T=5200 




log g =2.95 


log g =2.95 


log g = 2.50 


Mgi 


4.24 


4.39 


4.15 


Ah 


2.01 


2.15 


1.85 


Si I 


4.20 


4.38 


4.07 


Cai 


2.83 


2.98 


2.71 


Sen 


-0.29 


-0.20 


-0.54 


Til 


1.71 


1.83 


1.52 


Tin 


1.66 


1.74 


1.42 


Cri 


1.55 


1.68 


1.37 


Mm 


1.18 


1.29 


0.97 


Fei 


3.87 


4.00 


3.70 


Fen 


3.93 


4.02 


3.73 


Co I 


2.06 


2.20 


1.86 



Table C.2. Abundances for BS 16467-062 for different atmospheric parameters and the EWs of Lai et al. (2008). 



Ion 


ATLAS 


MARCS 




NOVER 








T=: 


5388 


T=: 


5200 




logg 


=3.04 


logg 


= 2.50 




A 


cr 


A 


cr 


Mgi 


4.06 


0.10 


3.98 


0.10 


Si I 


4.12 




4.01 




Cai 


2.90 


0.09 


2.78 


0.07 


Sen 


-0.30 


0.02 


-0.59 


0.02 


Til 


1.79 




1.57 




Tin 


1.71 


0.12 


1.43 


0.12 


Crn 


1.47 


0.06 


1.27 


0.07 


Mm 


1.05 


0.01 


0.82 


0.02 


Fei 


3.77 


0.19 


3.61 


0.18 


Fen 


3.76 


0.15 


3.54 


0.16 


Co I 


1.88 


0.06 


1.66 


0.06 


Nil 


2.80 


0.02 


2.60 


0.04 
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Table C.3. Abundances for BS 16467-062 from Paper V and the same model but EWs from Cohen et al. (2008), Lai et al. (2008). 



Ion 




EWs 






EWs 






Paper V 






Cohen et al. 2008 


Lai et al. 2008 












A 


cr 


N 


A 




N 


A 


a 


N 


Mgi 


3.99" 


0.11 


3 


3.98 


0.10 


4 


3.97 


0.09 


1 


Si I 


4.07 




1 


4.01 




1 


4.20 




1 


Cai 


2.71 


0.12 


3 


2.78 


0.07 


4 


2.94 


0.19 


12 


Sen 


-0.54 


0.03 


3 


-0.59 


0.02 


2 


-0.59 


0.06 


4 


Til 


1.52 


0.05 


2 


1.57 




1 


1.65 


0.17 


11 


Tin 


1.42 


0.11 


14 


1.43 


0.12 


9 


1.43 


0.18 


23 


Crn 


1.37 


0.10 


5 


1.27 


0.07 


4 


1.49 


0.29 


5 


Mm 


0.97 


0.22 


5 


0.82 


0.02 


2 


1.07 


0.03 


3 


Fei 


3.70 


0.16 


55 


3.61 


0.18 


52 


3.67 


0.13 


130 


Fe II 


3.73 


0.14 


8 


3.54 


0.16 


3 


3.79 


0.12 


4 


Co I 


1.86 


0.26 


3 


1.66 


0.06 


3 


1.70 


0.10 


4 


Nil 









2.60 


0.04 


2 


2.56 


0.03 


3 



" line 416.7 imi has been removed to compute the average 



